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Abstract 

A Chaitin Omega number is the halting probability of a universal Chaitin (self- 
delimiting Turing) machine. Every Omega number is both computably enumerable 
(the limit of a computable, increasing, converging sequence of rationals) and ran- 
dom (its binary expansion is an algorithmic random sequence). In particular, every 
Omega number is strongly non-computable. The aim of this paper is to describe a 
procedure, which combines Java programming and mathematical proofs, for com- 
puting the exact values of the first 64 bits of a Chaitin Omega: 

0000001000000100000110001000011010001111110010111011101000010000. 

Full description of programs and proofs will be given elsewhere. 



1 Introduction 

Any attempt to compute the uncomputable or to decide the undecidable is without doubt 



challenging, but hardly new (see, for example, Marxen and Buntrock |2^, Stewart |^ 
Casti |]TT|). This paper describes a hybrid procedure (which combines Java program- 
ming and mathematical proofs) for computing the exact values of the first 64 bits of a 
concrete Chaitin Omega number, Qu, the halting probability of the universal Chaitin 
(self-delimiting Turing) machine U, see |]13[. Note that any Omega number is not only 



uncomputable, but random, making the computing task even more demanding. 

Computing lower bounds for Qjj is not difficult: we just generate more and more 
halting programs. Are the bits produced by such a procedure exact? Hardly. If the 
first bit of the approximation happens to be 1, then sure, it is exact. However, if the 
provisional bit given by an approximation is 0, then, due to possible overflows, nothing 
prevents the first bit of Qu to be either or 1. This situation extends to other bits as 
well. Only an initial run of I's may give exact values for some bits of Qu- 

The paper is structured as follows. Section 2 introduces the basic notation. Com- 
putably enumerable (c.e.) reals, random reals and c.e. random reals are presented in 
Section 3. Various theoretical difficulties preventing the exact computation of any bits 
of an Omega number are discussed in Section 4. The register machine model of Chaitin 



T3| is discussed in Section 5. In section 6 we summarize our computational results con- 



cerning the halting programs of up to 84 bits long for U. They give a lower bound for 
flu which is proved to provide the exact values of the first 64 digits of Qu in Section 7. 

Chaitin has pointed out that the self-delimiting Turing machine constructed in 
the preliminary version of this paper ^ is universal in the sense of Turing (i.e., it is 
capable to simulate any self-delimiting Turing machine), but it is not universal in the 
sense of algorithmic information theory because the "price" of simulation is not bounded 
by an additive constant; hence, the halting probability is not an Omega number (but a 
c.e. real with some properties close to randomness). The construction presented in this 
paper is a self-delimiting Turing machine. Full details will appear in 



2 Notation 

We will use notation that is standard in algorithmic information theory; we will assume 
familiarity with Turing machine computations, computable and computably enumerable 
(c.e.) sets (see, for example. Bridges 0, Odifreddi |]2B|, Soare ||2^, Weihrauch [^) and 



elementary algorithmic information theory (see, for example, Calude 

By N, Q we denote the set of nonnegative integers (natural numbers) and rationals, 
respectively. If S* is a finite set, then #5* denotes the number of elements of S. Let 
S = {0, 1} denote the binary alphabet. Let S* be the set of (finite) binary strings, and 

the set of infinite binary sequences. The length of a string x is denoted by A 
subset A of S* is prefix-free if whenever s and t are in A and s is a prefix of t, then 
s = t. 

For a sequence x = xqXi ■ ■ - Xn - ■ ■ € S'^ and an nonnegative integer n > 1, x(?t,) 
denotes the initial segment of length n of x and Xi denotes the ith digit of x, i.e. x(?t,) = 
XqXi ■ ■ -Xn^i G S*. Due to Kraft's inequality, for every prefix-free set A C S*, Qa = 
X]se/i2~'''' lies in the interval [0,1]. In fact ^Ia is a probability: Pick, at random using 
the Lebesgue measure on [0, 1], a real a in the unit interval and note that the probability 
that some initial prefix of the binary expansion of a lies in the prefix-free set A is exactly 



Following Solovay |^D|, ^ we say that C is a ( Chaitin ) (self-delimiting Turing) ma- 



chine, shortly, a machine, if C is a Turing machine processing binary strings such that 
its program set (domain) PROGq = {x G S* | C{x) halts} is a prefix-free set of strings. 
Clearly, PROGc is c.e.; conversely, every prefix-free c.e. set of strings is the domain 
of some machine. The program-size complexity of the string x G S* (relatively to C) is 
Hc{x) = mm{\y\ \ y G E*, C{y) = x}, where min0 = oo. A major result of algorithmic 
information theory is the following invariance relation: we can effectively construct a 
machine U (called universal) such that for every machine C, there is a constant c > 
(depending upon U and C) such that for every x,y E S* with C{x) = y, there exists 
a string x' G E* with U{x') = y {U simulates C) and \x'\ < |x| -|- c (the overhead 
for simulation is no larger than an additive constant). In complexity-theoretic terms, 
Hu{x) < Hc{x) + c. Note that PROGjj is c.e. but not computable. 
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If C is a machine, then ^Ic = ^proGc represents its halting probabihty. When 
C = U is a universal machine, then its halting probability Qu is called a Chaitin 
number, shortly, Vt number. 



3 Computably Enumerable and Random Reals 

Reals will be written in binary, so we start by looking at random binary sequences. Two 
complexity-theoretic definitions can be used to define random sequences (see Chaitin 
IIT2I , [TtH): an infinite sequence x is Chaitin random if there is a constant c such that 
i^(x(n)) > n — c, for every integer n > 0, or, equivalently, lim„^oo -f^(x(n)) — n = 00. 
Other equivalent definitions include Martin-Lof [|^, ^ definition using statistical tests 
{Martin-Lof random sequences), Solovay measure-theoretic definition {Solovay ran- 



dom sequences) and Hertling and Weihrauch ||2T| topological approach to define random- 
ness [Hertling-Weihrauch random sequences). Independent proofs of the equivalence 
between Martin-Lof and Chaitin definitions have been obtained by Schnorr and Solovay, 



cf. [pTSj |I9|- In what follows we will simply call "random" a sequence satisfying one of 
the above equivalent conditions. Their equivalence motivates the following "randomness 
hypothesis" (Calude [^): A sequence is "algorithmically random" if it satisfies one of the 
above equivalent conditions. Of course, randomness implies strong non-computability 
(cf., for example, Calude but the converse is false. 

A real a is random if its binary expansion yi (i.e. a = O.xj is random. The choice 
of the binary base does not play any role, cf. Calude and Jiirgensen [|^, Hertling and 
Weihrauch ETI], Staiger |33]: randomness is a property of reals not of names of reals. 



Following Soare p8|, a real a is called c.e. if there is a computable, increasing 
sequence of rationals which converges {not necessarily computably) to a. We will start 
with several characterizations of c.e. reals (cf. Calude, Hertling, Khoussainov and Wang 
0). If O.y is the binary expansion of a real a with infinitely many ones, then a = 
E„ex. 2-"-i, where = \ = 1}. 

Theorem 1 Let a be a real in (0, 1]. The following conditions are equivalent: 

1. There is a computable, nondecreasing sequence of rationals which converges to a. 

2. The set {p G Q | p < a} of rationals less than a is c.e. 

3. There is an infinite prefix-free c.e. set ACS* with a = Qa- 

4. There is an infinite prefix-free computable set ACS* with a = Qa- 

5. There is a total computable function f : — > {0, 1} such that 

(a) If for some k, n we have f{k, n) = 1 and f{k, n-\-l) = then there is an I < k 
with f{l,n) = and f{l,n + 1) = 1 . 

(b) We have: k G <^=^ lim„^oo /(^, ''^) = 1- 
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We note that following Theorem |T], 5), given a computable approximation of a c.e. 
real a via a total computable function /, G X^^ <^=^ lim„^oo f{k,n) = 1; the values 
of f{k, n) may oscillate from to 1 and back; we will not be sure that they stabilized 
until 2^ changes have occurred (of course, there need not be so many changes, but in 
this case there is no guarantee of the exactness of the value of the /cth bit). 

Chaitin |]12| proved the following important result: 



Theorem 2 If U is a universal machine, then Vtu is c.e. and random. 

The converse of Theorem ^ is also true: it has been proved by Slaman based 
on work reported in Calude, Hertling, Khoussainov and Wang (see also Calude and 
Chaitin 0, Calude §, Downey 0): 



Theorem 3 Let a G (0, 1). The following conditions are equivalent: 

1. The real a is c.e. and random. 

2. For some universal machine U, a = Qu- 



4 The First Bits of An Omega Number 

We start by noting that 

Theorem 4 Given the first n bits of Qu one can decide whether U{x) halts or not on 
an arbitrary string x of length at most n. 

The first 10,000 bits of Vtjj include a tremendous amount of mathematical knowledge. 
In Bennett's words ]I|: 

[Vt] embodies an enormous amount of wisdom in a very small space 
. . . inasmuch as its first few thousands digits, which could be written on a 
small piece of paper, contain the answers to more mathematical questions 
than could be written down in the entire universe. 

Throughout history mystics and philosophers have sought a compact key to 
universal wisdom, a finite formula or text which, when known and understood, 
would provide the answer to every question. The use of the Bible, the Koran 
and the I Ching for divination and the tradition of the secret books of Hermes 
Trismegistus , and the medieval Jewish Cabala exemplify this belief or hope. 
Such sources of universal wisdom are traditionally protected from casual use 
by being hard to find, hard to understand when found, and dangerous to use, 
tending to answer more questions and deeper ones than the searcher wishes to 
ask. The esoteric book is, like God, simple yet undescribable. It is omniscient, 
and transforms all who know it . . . Omega is in many senses a cabalistic 
number. It can be known of, but not known, through human reason. To know 
it in detail, one would have to accept its uncomputable digit sequence on faith, 
like words of a sacred text. 
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It is worth noting that even if we get, by some kind of miracle, the first 10,000 digits 
of Qu, the task of solving the problems whose answers are embodied in these bits is 
computable but unrealistically difficult: the time it takes to find all halting programs of 
length less than n from O.Q0Q2 ■ ■ ■ ^n-i grows faster than any computable function of n. 

Computing some initial bits of an Omega number is even more difficult. According to 
Theorem |^, c.e. random reals can be coded by universal machines through their halting 
probabilities. How "good" or "bad" are these names? In ||12[| (see also [0, |1^), Chaitin 
proved the following: 

Theorem 5 Assume that ZFC^ is arithmetically sound^ Then, for every universal 
machine U, ZFC can determine the value of only finitely many bits of VLu- 

In fact one can give a bound on the number of bits of Vtjj which ZFC can determine; 



this bound can be explicitly formulated, but it is not computable. For example, in |T5 
Chaitin described, in a dialect of Lisp, a universal machine U and a theory T, and 
proved that U can determine the value of at most H{T) + 15, 328 bits of Qu] H{T) is 
the program-size complexity of the theory T, an uncomputable number. 

Fix a universal machine U and consider all statements of the form 

"The ra*^ binary digit of the expansion ofilu is fc", (1) 

for all n > 0,k = 0,1. How many theorems of the form (|l]) can ZFC prove? More 
precisely, is there a bound on the set of non-negative integers n such that ZFC proves a 
theorem of the form (|l])? From Theorem ^ we deduce that ZFC can prove only finitely 
many (true) statements of the form (p. This is Chaitin information-theoretic version of 



Godel's incompleteness (see [15, ^ 



Theorem 6 // ZFC is arithmetically sound and U is a universal machine, then almost 
all true statements of the form are unprovable in ZFC. 

Again, a bound can be explicitly found, but not effectively computed. Of course, for 
every c.e. random real a we can construct a universal machine U such that a = Qu and 
ZFC is able to determine finitely (but as many as we want) bits of Qu. 

A machine U for which Peano Arithmetic can prove its universality and ZFC cannot 
determine more than the initial block of 1 bits of the binary expansion of its halting 
probability, Qu, will be called Solovay machine^ To make things worse Calude 
proved the following result: 

Theorem 7 Assume that ZFC is arithmetically sound. Then, every c.e. random real 
is the halting probability of a Solovay machine. 



^Zermelo set theory with choice. 

^That is, any theorem of arithmetic proved by ZFC is true. 
3 Clearly, U depends on ZFC. 
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For example, if a G (3/4, 7/8) is c.e. and random, then in the worst case ZFC can 
determine its first two bits (11), but no more. For a G (0, 1/2) we obtained Solovay's 
Theorem [pl[] : 



Theorem 8 Assume that ZFC is arithmetically sound. Then, every c.e. random real 
a G (0, 1/2) is the halting probability of a Solovay machine which cannot determine any 
single bit of a. No c.e. random real a G (1/2, 1) has the above property. 

The conclusion is that the worst fears discussed in the first section proved to mate- 
rialize: In general only the initial run of I's (if any) can be exactly computed. 



5 Register Machine Programs 



First we start with the register machine model used by Chaitin |T^. Recall that any 
register machine has a finite number of registers, each of which may contain an arbitrarily 
large non-negative integer. The list of instructions is given below in two forms: our 



compact form and its corresponding Chaitin |]T3[ version. The main difference between 
Chaitin's implementation and ours is in the encoding: we use 7 bit codes instead of 8 
bit codes. 

L: ? LI (L: GOTO LI) 

This is an unconditional branch to LI. LI is a label of some instruction in the program 
of the register machine. 

L: A R LI (L: JUMP R LI) 

Set the register R to be the label of the next instruction and go to the instruction 
with label LI. 

L: @ R (L: GOBACK R) 

Go to the instruction with a label which is in R. This instruction will be used in 
conjunction with the jump instruction to return from a subroutine. The instruction 
is illegal (i.e., run-time error occurs) if R has not been explicitly set to a valid label 
of an instruction in the program. 

L: = Rl R2 LI (L: EQ Rl R2 LI) 

This is a conditional branch. The last 7 bits of register Rl are compared with the 
last 7 bits of register R2. If they are equal, then the execution continues at the 
instruction with label LI. If they are not equal, then execution continues with the 
next instruction in sequential order. R2 may be replaced by a constant which can be 
represented by a 7-bit ASCII code, i.e. a constant from to 127. 
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L: # Rl R2 LI 



(L: NEQ Rl R2 LI) 



This is a conditional branch. The last 7 bits of register Rl are compared with the 
last 7 bits of register R2. If they are not equal, then the execution continues at 
the instruction with label LI. If they are equal, then execution continues with the 
next instruction in sequential order. R2 may be replaced by a constant which can be 
represented by a 7-bit ASCII code, i.e. a constant from to 127. 

L: ) R (L: RIGHT R) 



Shift register R right 7 bits, i.e., the last character in R is deleted. 
L: ( Rl R2 (L: LEFT Rl R2) 

Shift register Rl left 7 bits, add to it the rightmost 7 bits of register R2, and then 
shift register R2 right 7 bits. The register R2 may be replaced by a constant from 
to 127. 

L: & Rl R2 (L: SET Rl R2) 

The contents of register Rl is replaced by the contents of register R2. R2 may be 
replaced by a constant from to 127. 

L: ! R (L: READ R) 

One bit is read into the register R, so the numerical value of R becomes either or 
1. Any attempt to read past the last data- bit results in a run-time error. 

L: / (L: DUMP) 

All register names and their contents, as bit strings, are written out. This instruction 
is also used for debugging. 

L: % (L: HALT) 



Halts the execution. This is the last instruction for each register machine program. 



A register machine program consists of finite list of labeled instructions from the 
above list, with the restriction that the HALT instruction appears only once, as the last 
instruction of the list. The data (a binary string) follows immediately the HALT instruc- 
tion. The use of undefined variables is a run-time error. A program not reading the 
whole data or attempting to read past the last data-bit results in a run-time error. Be- 
cause of the position of the HALT instruction and the specific way data is read, register 
machine programs are Chaitin machines. 
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To be more precise, we present a context-free grammar G — {N, E, P, S) in Backus- 
Naur form which generates the register machine programs. 



1^1) N is the finite set of nonterminal variables: 

N = {SjUlNSTUTOKEN 

INST = {(RMSi^,),(?i^g),(Ain3),(@i^g),(=In3),(#I^g), 
Oins)' (dns)' (-^Ins)' Clns)' (/ins)' (%ns)} 

TOKEN = {(Data), (Label), (Register), (Constant), 
(Special), (Space), (Alpha), (LS)} 



(2) E, the alphabet of the register machine programs, is a finite set of terminals, disjoint 
from N: 

S = (Alpha) u (Special) u (Space) u (Constant) 

(Alpha) = {a,b,c,..., z} 

(Special) = {:,/,?, A, ©,=,#,), (,&,!,?,%} 

(Space) = {'space', 'tab'} 

(Constant) = {d\0<d< 127} 



(3) P (a subset of x (NUT,)*) is the finite set of rules (productions): 

S ^ (RMSins)*(%i^g)(DATA) 
(0|1)* 

I (1|2| ... |9)(0|1|2| ... |9)* 
: (Space)* 

(Alpha) ((Alpha) u (0|1|2| . . . |9))* 



(Data) 
(Label) 
(LS) 
(Register) 
(RMSins) 



(•Ins) I i^lns) I (®Ins) I (=Ins) I (#Ins) 
Oins) I (dns) I (%ns) I Clns) I (/ins) 



(L: HALT) 
(%Ins) ^ (Label) (LS)% 

(L: GOTO LI) 

(•Ins) ^ (Label) (LS)?(Space)* (Label) 

(L: JUMP R LI) 

(Aj^g) (Label)(LS) A (Space)*(Register)(Space)+(Label) 

(L: GOBACK R) 

(®Ins) ^ (Label) (LS)@(Space)* (Register) 



(L: EQ R 0/127 LI or L: EQ R R2 LI) 

(=Ins) ~^ (Label)(LS) = (Space)*(Register)(Space)+(Constant)(Space)"^ 
(Label) | (Label) (LS) = (Space)* (Register) (Space)^ (Register) 
(Space)^ (Label) 

(L: NEQ R 0/127 LI or L: NEQ R R2 LI) 

(#jj^g) (Label) (LS)#(Space)* (Register) (Space)^ (Constant) (Space)'^ 
(Label) | (Label) (LS)#(Spage)* (Register) (Space) + 
(Register) (Space)^ (Label) 



(L: RIGHT R) 

Oins) ~^ (Label) (LS)) (Space)* (Register) 

(L: LEFT R LI) 

((ins) ~^ (Label)(LS)((Space)*(Register)(Space)^(Constant) 
(Label) (LS) ((Space)* (Register) (Spage)"^ (Register) 

(L: SET R 0/127 or L: SET R R2) 

(&j^g) (Label) (LS)&(Space)* (Register) (Space)^ (Constant) 
(Label) (LS)&;(Space)* (Register) (Space)"^ (Register) 

(L: READ R) 

Clns) ~^ (Label) (LS)! (Space)* (Register) 

(L: DUMP) 

(/ins) ^ (Label)(LS)/ 



(4) S E N is the start symbol for the set of register machine programs. 

It is important to observe that the above construction is universal in the sense of 
algorithmic information theory (see the discussion of the end of Section 1). Register 
machine programs are self-delimiting because the HALT instruction is at the end of any 
valid program. Note that the data, which immediately follows the HALT instruction, is 
read bit by bit with no endmarker. This type of construction has been first programmed 
in Lisp by Chaitin |]T3|, IT? . 

To minimize the number of programs of a given length that need to be simulated, 
we have used "canonical programs" instead of general register machines programs. 
A canonical program is a register machine program in which (1) labels appear in 
increasing numerical order starting with 0, (2) new register names appear in increasing 
lexicographical order starting from 'a', (3) there are no leading or trailing spaces, (4) 
operands are separated by a single space, (5) there is no space after labels or operators, 
(6) instructions are separated by a single space. Note that for every register machine 
program there is a unique canonical program which is equivalent to it, that is, both 
programs have the same domain and produce the same output on a given input. If x is 
a program and y is its canonical program, then \y\ < \x\. 
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Here is an example of a canonical program: 



0:!a l:"b 4 2:!c 3:?11 4:=a 8 5:&c 110 6:(c 101 7:@b 8:&c 101 
9: (c 113 10:0b 11:7.10 



To facilitate the understanding of the code we rewrite the instructions with additional 
comments and spaces: 

: ! a // read the first data bit into register a 

1 : " b 4 // jump to a subroutine at line 4 

2 : ! c //oil return from the subroutine call c is written out 

3 : ? 11 // go to the halting instruction 

4 : = a 8 // the right most 7 bits are compared with 127; if they 

/ / are equal, then go to label 8 

5:& c 'n' // else, continue here and 

6: ( c 'e' // store the character string 'ne' in register c 

7 : @ b //go back to the instruction with label 2 stored in 

/ / register b 

8:& c 'e' // store the character string 'eq' in register c 
9: ( c 'q' 
10:@ b 

11 :°/o // the halting instruction 

10 // the input data 

For optimization reasons, our particular implementation designates the first maximal 
sequence of SET/LET instructions as (static) register pre-loading instructions. We "com- 
press" these canonical programs by (a) deleting all labels, spaces and the colon symbol 
with the first non-static instruction having an implicit label 0, (b) separating multiple 
operands by a single comma symbol, (c) replacing constants with their ASCII numerical 
values. The compressed format of the above program is 



!a"b,4!c?ll=a,0,8&c,110(,c,101(§b&,c,101(,c,113@by.l0 



Note that compressed programs are canonical programs because during the process 
of "compression" everything remains the same except for the elimination of space. Com- 
pressed programs use an alphabet with 49 symbols (including the halting character). The 
length is calculated as the sum of the program length and the data length (7 times the 
number of characters). For example, the length of the above program is 7 x (49-1-2) = 357. 

For the remainder of this paper we will be focusing on compressed programs. 



6 Solving the Halting Problem for Programs Up to 
84 Bits 

A Java version interpreter for register machine compressed programs has been imple- 
mented; it imitates Chaitin's universal machine in |jl3|]. This interpreter has been used 
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to test the Halting Problem for all register machine programs of at most 84 bits long. 
The results have been obtained according to the following procedure: 

1. Start by generating all programs of 7 bits and test which of them stops. All strings 
of length 7 which can be extended to programs are considered prefixes for possible 
halting programs of length 14 or longer; they will simply be called prefixes. In 
general, all strings of length n which can be extended to programs are prefixes 
for possible halting programs of length n + 7 or longer. Compressed prefixes are 
prefixes of compressed (canonical) programs. 

2. Testing the Halting Problem for programs of length n G {7, 14, 21, . . . , 84} was 
done by running all candidates (that is, programs of length n which are exten- 
sions of prefixes of length n — 7) for up to 100 instructions, and proving that 
any generated program which does not halt after running 100 instructions never 
halts. For example, (uncompressed) programs that match the regular expression 
"0:\" a 5.* 5:\? 0" never halt on any input. 

For example, the programs " lalblalb/ZlOlOlOlO" and " ! a?oyol0101010" produce 
a run-time errors; the first program 'under reads' the data and the second one 'over 
reads' the data. The program " ! a?l IbZlOlO" loops. 

One would naturally want to know the shortest program that halts with more than 
100 steps. If this program is larger than 84 bits, then all of our looping programs never 
halt. The trivial program with a sequence of 100 dump instructions runs for 101 steps 
but can we do better? The answer is yes. The following family of programs {Pi, P2, . . .} 
recursively count to 2* but have linear growth in size. The programs Pi through P4 are 
given below:0 

/&a,0=a,l,5&a,l?2y„ 
/&a,0&b,0=b,l,6&b,l?3=a,l,9&a,l?2°/o 

/&a , O&b , O&c , 0=c , 1 , 7&C , l?4=b , 1 , lO&b , 1 ?3=a , 1 , 13&a , l?2l 

/&a , O&b , O&c , O&d , 0=d , 1 , 8&d , l?5=c , 1 , 1 l&c , l?4=b , 1 , 14&b , 1 ?3=a , 1 , 17&a , 172% 

In order to create the program Pj+i from Pj only 4 instructions are added, while 
updating 'goto' labels. 

The running time t{i), excluding the halt instruction, of program Pj is found by the 
following recurrence: t{l) = 6, t{i) = 2-t(i — 1) + 4. Thus, since t (4) = 86 andt(5) = 156, 
P5 is the smallest program in this family to exceed 100 steps. The size of P5 is 86 bytes 
(602 bits), which is smaller than the trivial dump program of 707 bits. It is an open 
question on what is the smallest program that halts after 100 steps. An hybrid program, 
given below, created by combining P2 and the trivial dump programs is the smallest 
known. 

&a, 0/&b , 0/////////////////////=b , 1 , 26&b , l?2=a, 1 , 29&a , l?Ol 

Hn all cases the data length is zero. 
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This program of 57 bytes (399 bits) runs for 102 steps. 

Note that the problem of finding the smallest program with the above property is 



undecidable (see 16 



The distribution of halting compressed programs of up to 84 bits for U, the universal 
machine processing compressed programs, is presented in Table 1. All binary strings 
representing programs have the length divisible by 7. 



Program plus 


Number of 


Program plus 


Number of 


data length 


halting programs 


data length 


halting programs 


7 


1 


49 


1012 


14 


1 


56 


4382 


21 


3 


63 


19164 


28 


8 


70 


99785 


35 


50 


77 


515279 


42 


311 


84 


2559837 



Table 1. Distribution of halting programs 



7 The First 64 Bits of Qu 

Computing all halting programs of up to 84 bits for U seems to give the exact values of 
the first 84 bits of Qu. False! To understand the point let's first ask ourselves whether 
the converse implication in Theorem ^ is true? The answer is negative. Globally, if we 
can compute all bits of Qu, then we can decide the Halting Problem for every program for 
U and conversely. However, if we can solve for U the Halting Problem for all programs 
up to bits long we might not still get any exact value for any bit of Qu (less all values 
for the first bits). Reason: A large set of very long halting programs can contribute 
to the values of more significant bits of the expansion of Qu- 

So, to be able to compute the exact values of the first A^ bits of Qu we need to be 
able to prove that longer programs do not affect the first A^ bits of Qu- And, fortunately, 
this is the case for our computation. Due to our specific procedure for solving the 
Halting Problem discussed in Section 6, any compressed halting program of length n 
has a compressed prefix of length n — 7. This gives an upper bound for the number of 
possible compressed halting programs of length n. 

Let ^l^ be the approximation of Qu given by the summation of all halting programs 
of up to n bits in length. 

Compressed prefixes are partitioned into two cases — ones with an HALT (%) instruc- 
tion and ones without. Hence, halting programs may have one of the following two forms: 
either "x y HALT u" , where a; is a prefix of length k not containing HALT, y is a sequence 
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of instructions of length n — k not containing HALT and u is the data of length m > 0, or 
"x u" , where x is a prefix of length k containing one occurrence of HALT followed by data 
(possibly empty) and u is the data of length m > 1. In both cases the prefix x has been 
extended by at least one character. Accordingly, the "tail" contribution to the value of 

oo 

^u = T. E 2-1-1 

n=0{[^„|=„_[/(^) halts} 

is bounded from above by the sum of the following two convergent series (which reduce 
to two independent sums of geometric progressions): 

oo oo 

J2 J2 I prefix X not containing HALT, \x\ = k} ■ 47""'= • 1 ■2" .l28"("+™+^), 

m=On=k'' ^ ' y jjALT « 

and 

oo 

#{x I prefix x containing HALT, |x| = k} ■ 2"^ ■128'^"'^''\ 

The number 47 comes from the fact that the alphabet has 48 characters and the last 
instruction before the data is HALT (%). 

There are 402906842 prefixes not containing HALT and 1748380 prefixes containing 
HALT. Hence, the "tail" contribution of all programs of length 91 or greater is bounded 
by: 

oo oo oo 

-(m+13) 



m=0 n=13 m=0 



J2 E 402906842 ■ 47""^^ ■ 2"^ ■ i28~("+™+i) + ^ 1748380 ■ 2"^ ■ 128" 

m=0 

°° / 47 \ " 1 

y {—) + 1748380 < 2-^^ (2) 

V128/ 63- 12813 ' ^ ^ 



64 

402906842 



that is, the first 68 bits of Qfj "may be" correct by our method. Actually we do not have 
68 correct bits, but only 64 because adding a 1 to the 68th bit may cause an overflow up 
to the 65th bit. From (H) it follows that no other overflows may occur. 



The following list presents the main results of the computation: 

= 0.0000001 
njf = 0.00000010000001 

f]2i = 0.000000100000010000011 

nff = 0.0000001000000100000110001000 

nff = 0.00000010000001000001100010000110010 

nf^ = 0.000000100000010000011000100001101000110111 

nff = 0.0000001000000100000110001000011010001111101110100 

f]56 ^ 0.00000010000001000001100010000110100011111100101100011110 

f]63 ^ 0.000000100000010000011000100001101000111111001011101100111011100 
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fi7o ^Q QQQQQQ^QQQQQQ^OOOOOllOOOlOOOOllOlOOOllllllOOlOlllOlllOOlllOOllllOOlOOl 
= 0.0000001000000100000110001000011010001111110010111011101000001110000010 
1001111 

f]84^ 0.0000001000000100000110001000011010001111110010111011101000010000011110 
11011011011101 



The exact bits are underlined in the 84 approximation: 

= 0. 0000001000000100000110001000011010001111110010111011101000010000 011110 
11011011011101 



In summary, the first 64 exact bits of ^lu are: 

0000001000000100000110001000011010001111110010111011101000010000 



8 Conclusions 

The computation described in this paper is the first attempt to compute some initial 
exacts bits of a random real. The method, which combines programming with mathe- 
matical proofs, can be improved in many respects. However, due to the impossibility of 
testing that long looping programs never actually halt (the undecidability of the Halting 
Problem), the method is essentially non-scalable. 

As we have already mentioned, solving the Halting Problem for programs of up to n 
bits might not be enough to compute exactly the first n bits of the halting probability. 
In our case, we have solved the Halting Problem for programs of at most 84 bits, but we 
have obtained only 64 exact initial bits of the halting probability. 

Finally, there is no contradiction between Theorem ^ and the main result of this 
paper. Jl's are halting probabilities of Chaitin universal machines, and each Q is the 
halting probability of an infinite number of such machines. Among them, there are those 
(called Solovay machines in [^) which are in a sense "bad" as ZFC cannot determine 
more than the initial run of I's of their halting probabilities. But the same VL can be 
defined as the halting probability of a Chaitin universal machine which is not a Solovay 
machine, so ZFC, if supplied with that different machine, may be able to compute more 
(but always, as Chaitin proved, only finitely many) digits of the same VL. Such a machine 
has been used for the VL discussed in this paper. 



The web site [ftp : //f tp . cs . auckland. ac .nz/pub/CDMTCS/Omega/| contains all pro- 



grams used for the computation as well as all intermediate and final data files (3 giga- 
bytes in gzip format). 
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